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The  conventional  approach  to  ship  steering  is  to  regard  the  ship  as  a  single  input, 
single  output  system  without  cross-coupling  or  interaction  between  speed,  yaw  and  roll. 
This  approach  has  found  successful  application,  particularly  in  conventional  vessels 
where  the  amount  of  cross-coupling  is  normally  slight.  But,  as  a  result  of  tight 
maneuvering,  the  modem  warship  suffers  severe  cross-coupiing  effects  because  of  large 
control  surfaces,  high  speed  and  low  tonnage.  Consequently,  the  adoption  of  a 
multivariable  approach  to  ship  steering  would  appear  to  be  more  suited  for  the  design 
of  a  steering  control  system. 

This  thesis  describes  the  results  of  a  simulation  study  of  pre-compensator  design 
to  suppress  the  undesirable  cross-coupling  effects  between  speed,  yaw  and  roil. 

Simulation  studies  using  DSL  and  Function  Minimization  are  the  basis  for 
accomplishing  the  design. 

Simulation  results  presented  indicate  that  the  adoption  of  a  multi-input,  multi¬ 
output  approach  would  result  in  a  significant  improvement  in  the  combined  steering 
and  stabilization  problem  of  a  warship.  j7 .•  ,  . .  ,  •  -  ' 
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I.  INTRODUCTION 


Modern  warships  must  be  highly  maneuverable  to  satisfy  numerous  operational 
requirements.  However,  the  steering  characteristics  of  a  modern  ship  are  nonlinear,  so 
severe  interaction  or  cross-coupling  exists  between  the  control  surface  inputs  and  the 
controlled  outputs.  Therefore,  warship  steering  is  a  complex  multivariable  control 
problem. 

Interaction  between  roll,  yaw  and  speed  are  pronounced  in  warships  because  of 
their  length  to  beam  ratio  and  relatively  large  control  surfaces.  However,  warships 
should  be  able  to  execute  high  speed  maneuvers  while  maintaining  their  fightina 
capability.  This  is  generally  not  possible  due  to  severe  cross-coupiing. 

In  this  thesis,  the  development  ot  a  pre-compensator  to  reduce  the  cross-coupiing 
effects  which  are  present  :n  the  steering  characteristics  of  a  modern  ship  is  introduced. 
The  design  method  adopted  by  Roberts  [Ref.  lj  uses  the  Direct  Nvquist  Array  (DMA) 
frequency  response  technique  as  defined  by  Rosenbrock  [Ref.  2|  and  Fricker  [Ref.  3]. 
This  thesis  uses  the  basic  ship  model  and  the  results  which  were  reported  bv  Roberts 
which  have  the  potential  to  produce  improved  seakeeping  and  ship  stability. 
Emphasis  in  this  thesis  is  placed  on  optimizing  the  parameters  of  a  pre-compensator 
using  Function  Minimization  (F.M.)  via  a  digital  computer. 

Simulation  studies  employed  Function  Minimization  techniques  together  with  the 
Dynamic  Simulation  Language  (DSL")  package. 
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II.  BASIC  MODEL 


The  model  and  data  used  in  the  study  is  that  proposed  by  Roberts  [Ref.  I).  The 
structure  of  the  ship  model  is  shown  in  Figure  2.1  and  the  elements  of  the  transfer 
function  matrix,  G(S),  which  were  derived  using  curve  fitting  techniques  to  measured 
step  and  frequency  response  data,  are  given  in  Table  1. 


Figure  2.1  Multivariable  Structure  of  a  Warship  Model. 


The  nonlinear  nature  of  the  ship  dynamics  is  demonstrated  by  the  change  in 
steady  state  gain  parameters  as  shown  in  Table  2. 
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For  ease  of  anahsis  the  system  inputs  and  outputs  are  used  as  defined  below  in 
Table  3. 


TABLE  3 

MAX  VALUES  OF  INPUT  AND  OUTPUT 


-  '  . . .  -  —  —  ■  - 

INPUTS 

MAX  VALUE 

Fin  Angle 

±27* 

Rudder  Angle 

±30* 

Power 

±10% 

OUTPUTS 

MAX  VALUE 

Roll  Angle 

±15* 

Yaw  Angle 

±120* 

Forward  Speed 

30  Kts 

Figure  2.2  and  Figure  2.3  show  the  time  responses  for  the  ship  model  when  we 
use  10%  demands  for  rudder  and  fin  at  12  kts. 

As  can  be  seen  in  Figure  2.2,  there  are  pronounced  cross-coupling  effects  in  both 
roll  and  speed  outputs. 

In  this  thesis,  10%  demand  for  rudder  at  12  kts  will  be  always  used  for 
convenience. 
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III.  PRE-COMPENSATOR  DESIGN 


The  aim  of  including  a  pre-compensator,  K(S),  is  to  decouple  the  interaction 
present  in  G(S),  thus  enabling  the  reduction  of  interaction  or  cross-coupling  in  the 
system.  The  action  of  the  pre-compensator  is  to  propagate  the  three  input  demands  in 
such  a  way  that  each  input  affects  its  associated  output  only. 

This  chapter  was  extracted  from  Roberts  [Ref.  1). 

The  system  configuration  with  the  pre-compensator  included  is  shown  in  Figure 
3.1. 


The  approach  used  in  this  study  is  to  develop  a  pre-compensator,  K(S),  which 
totally  diagonalizes  the  pair  G(S)K(S)  (  =  Q(S)  ).  This  is  in  effect  non-interacting 
control  and  can  result  in  the  elements  of  K(S)  having  high  order  which  may  prove  a 
problem  when  it  comes  to  implementation,  particularly  if  this  is  to  be  achieved  using 
analog  techniques.  If  necessary  the  complexity  of  the  elements  of  K(S)  can  be  reduced 
using  standard  reduction  routines  while  maintaining  diagonal  dominance  in  G(S)  K(S). 


The  diagonalizing  pre-compensator  was  produced  using  the  method  suggested  by 
Flicker  [Ref.  3].  This  technique  produces  an  initial  "ideal"  pre-compensator  K(S)  which 
diagonalizes  Q(S).  The  elements  of  K(S),  which  can  be  high  order,  are  then  reduced  to 
a  simpler  form  while  ensuring  that  diagonal  dominance  in  Q(S)  is  maintained.  This 
involves  expressing  G(S)  as: 


G(S)  -  R(S)  A(S) 


(eqn  3. 


where  R(S)  is  a  diagonal  matrix  formed  by  extracting  common  row  elements  from  G(S') 
so  that  A(S)  contains  numerator  polynomial  elements  only. 

The  decoupling  pre-compensator  is  therefore  given  by: 


K(S)  -  A(S)-1 


(eqn  3.2', 


However,  this  method  is  only  possible  if  |A(s)|  =  0  has  all  stable  factors,  i.e.,  all 
roots  of  the  characteristic  equation  are  in  the  left  half  plane.  If  this  is  not  the  case, 
then  K(S)  is  formed  from: 


K(S)  - 


Adjoint  A(S) 
Realization  Factor 


(eqn  3.3) 


where  the  realization  factor  will  contain  the  stable  factors  of  |A(S)|  together  with  other 
suitable  lag  elements  to  make  K(S)  physically  realizable. 

The  elements  of  the  diagonal  matrix  R(S)  and  the  elements  of  A(S)  so  formed  are 
given  in  Table  4  and  Table  5  respectively.  As  |A(S)|  contains  all  stable  factors,  the 
ideal  pre-compensator  can  be  formed  directly  from  Equation  3.2.  The  results  of  this 
operation  are  given  in  Table  6. 

It  is  necessary  to  scale  the  columns  of  K(S)  to  ensure  that  the  steady-state  gains 
of  the  diagonal  elements  of  Q(S)  remain  the  same  as  those  of  G(S).  The  elements  of 
this  final  pre-compensator  are  given  in  Table  7.  Table  8  gives  the  speed  related  gain 
variation. 

Figure  3.2  and  Figure  3.3  show  that  the  addition  of  the  pre-compensator  has 
resulted  in  an  improvement  in  outputs  at  12  kts. 


18 


Cvnow 


WWW 


nv^v  mvwj  w  >r 


i 


TABLE  4 

ELEMENTS  OF  DIAGONAL  MATRIX  R(S) 


|Rn(S) 

I 

'r22(S) 

R3jlS> 


1 

(  S2  +  .  06S  +  .  25  )  (  S2  +  .  02S  +  .  154)(8.2S  +  1) 

_ 1 _ 

S(  .  434S  +  1 )  (  6.  62S  -  1)(4.  1SS  +  1) 


(  S2  +  .  199S  -r  .  99  )  (  10S  +  1 )  (  24S  *  1) 


i 


TABLE  5 

ELEMENTS  OF  MATRIX  A(S) 


AU(S) 

= 

.25KU(S2  +  .  2S  + .  15  )  (  8.  2S  +  1) 

Ai2(s) 

s  - 

.15K.2(S2  +  .  06S  ♦  .  25  ) (  -8.  57S  +  1) 

a^s) 

= 

K22 

A31(  S) 

= 

.  1K31(  10S  +  1)(  10.  053  +  1) 

A.;(  3) 

= 

K32(  32  +  .  2S  ♦  .1) 

■  A33(  S) 

= 

K,3(  S2  ♦  .  2S  ♦  .  1)(  10.  05S  -  1 ) 

a;3;  3 ) 

= 

a21(  S )  =  A23(  S)  =  0 

|A( S) 1 

1 

s 

.  25K,.K,,K,,(  S2  +  .  2S  +.15)(S2  ♦  .  2S  ♦  .1) 
ll  33(  8.  23S  *  1)(  10.  05S  ♦  1) 

This  "ideal"  pre-compensator  can  be  reduced  to  the  individual  elements  of  Kisi 
using  step  response  data. 

After  the  reduction  process  is  completed,  the  elements  of  the  reduced  pre- 
compensator  so  formed  are  given  in  Table  9.  The  terms  Kcl,  Kc2>  and  Kc3  given  in 
Tabie  9  are  the  speed  related  compensator  gain  changes  necessary  for  the  pre- 
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TABLE  6 

ELEMENTS  OF  DIAGONALIZING  PRE-COMPENSATOR 


KU(S) 

!k,,<  s) 

1  k:2(  s  ) 
K31(S) 

;  k32(  s ) 

;:<33(S) 

,  k13(S) 


i 


_  _ .  00649 _  I 

53.  19S3  +  17.  136S2  +  9.  5S  +  1  | 

_  -.  00454(  -34.  143S3  +  2S2  -  8.  29S  +1)  j 

"  53.19S3  ~  17.  136S2  +  9.  5S  +  1  ! 

=  . 00285 

. 00376( 10S  +  1)  | 

"  534S5  +  278S4  +  182.  6S3  +  46.  15S2  ~  11.  5S  +  1  ! 

_  .  00013(  80050S5  ♦  15650S4  -  20665S3  +  2255S2  +  3S  +  1) 

5352S0  -  3326S5  +  2113S4  +  S46S3  +  162S2  *  21.  S3  *  1! 

. 00285 

100S3  +30S2  +12S  *1 

=  K;I(S)  =  K:3(S)  =  0  ' 


compensator  to  cope  with  the  non-linearity  of  the  warship  model  and  these  are  the 
same  as  in  Table  8. 

Figure  3.4  and  Figure  3.5  show  that  interactions  between  inputs  and  outputs 
have  not  oeen  eliminated  in  the  system's  time  response  when  we  use  the  reduced  order 
pre-compensator  aiso. 
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TABLE  7 

FINAL  ELEMENTS  OF  DIAGONALIZING  PRE-COMPENSATOR 


Kn(S) 


k12(S) 
K:2(  S ) 
:<31(S) 

:<3;(  s ) 

K33(S) 

k13(S) 


1 

53.  19S3  +  17.136S2  +  9.  5S  +  1 
K.[( -34.  143S3  +  2S2  -  8.  29S  +  1) 

53.  19S3  +  17.-136S2  +  9.  5S  +  1 
1 

_  Kc2(10S  f  1) _ 

534S5  +  278S4  +  182.  6S3  +  46.  15S2  11.  5S  +  1 

K„3(  80050S5  +  15650S4  +  20665S3  +  2255S2  +  3S  -  1) 
5352S0  +  3326S5  +  2113S4  +  646S3  +  162S2  +  21.  6S 
1 

100S3  +  30S2  +12S  +1 
K^S)  =  IC3(S)  =  0 


TABLE  8 

PRE-COMPENSATOR  GAIN  VARIATION 


TABLE  9 

FINAL  REDUCED  PRE-COMPENSATOR  ELEMENTS 


IV.  COMPUTER  SIMULATION 


From  the  results  discussed  so  far  (all  of  which  are  from  Roberts  but  were 
repeated  at  N.P.S.  as  part  of  this  thesis),  it  can  be  seen  that  the  pre-compensator 
substantially  reduces  the  interaction  intensity  of  the  compensated  system  and  results  in 
improved  ship  stability  and  minimal  loss  of  speed  while  executing  normal  maneuvers. 
Also,  it  has  been  shown  that  the  outputs  of  ship's  roil  angle,  heading  angle  and  speed 
can  be  controlled  by  a  low-order  pre-compensator.  However,  it  can  also  be  seen  that 
large  cross-coupling  effects  exist  which  cause  long  settling  times.  Therefore,  a  method 
must  still  be  devised  which  controls  the  interactions  between  inputs  and  outputs  and 
shortens  settling  time.  The  outputs  of  ship's  roll  angle,  heading  angie  and  speed  can  be 
changed  by  a  few  factors. 

First  of  ail.  this  thesis  will  introduce  the  philosophy  of  Function  Minimization 
(F.M.),  which  is  the  main  theme  of  this  paper  and  introduce  a  few  factors,  which  have 
important  effects  on  outputs. 

Next,  using  above  factors,  this  thesis  will  design  other  pre-compensators  and 
simulate  them. 

A.  PHILOSOPHY  OF  FUNCTION  MINIMIZATION 

Classical  control  theory  has  historically  been  applied  to  the  design  problem  with 
the  assistance  of  graphical  presentations  and  trial  and  error  methods.  Such  methods 
have  been  quite  successful  in  the  development  of  good  control  systems,  but  do  not 
answer  the  question  "Is  this  the  best  system  possible?". 

If  a  given  function  has  a  minimum  within  the  range  of  permitted  parameter 
variation,  there  exist  numerical  methods  which  can  be  used  to  find  the  minimum. 
These  numerical  methods  have  been  programmed  and  most  computer  libraries  contain 
one  or  more  subroutines  which  can  be  used  for  Function  Minimization.  This  thesis 
uses  the  HOOKE  [Ref.  4]  subroutine  of  DSL. 

Since  one  has  the  freedom  to  select  the  cost  function  to  be  used,  certainly  any  of 
those  used  within  optimal  control  theory  can  be  chosen.  We  can  therefore  design  an 
"optimal  controller"  without  using  the  conventional  theoretical  approach. 

We  can  also  use  any  of  the  well  known  performance  indices  as  a  cost  function, 
i.e.,  jE2  dt,  J|Ej  dt,  Jj Ejt  dt  are  easily  evaluated  and  minimized. 
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In  addition,  one  can  select  the  cost  function  to  suit  the  particular  specification  of 
the  problem. 

Figure  4.1  shows  the  block  diagram  for  simulation  using  Function  Minimization. 


Figure  4.1  Block  Diagram  of  Function  Minimization. 


B.  CHOICE  OF  COMPENSATORS 

Use  of  Function  Minimization  in  the  computer  is  relatively  expensive,  so  some 
preliminary  analysis  and  design  is  desirable  to  avoid  excessive  computer  time.  For  a 
simple  system  (  one  output  )  a  BODE  design  or  a  ROOT  LOCUS  design  might  be  a 
good  starting  point.  The  compensator  thus  found  is  then  simulated  and  improved 
through  use  of  Function  Minimization.  For  the  ship  problem  in  this  thesis,  the 
preliminary  work  was  done  by  Roberts  as  can  be  seen  in  Chapters  II  and  III. 
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This  thesis  uses  Roberts'  compensators  and  outputs  as  a  starting  point.  The 
desired  output  curves  were  obtained  by  modifying  the  known  outputs.  The  questions 
to  be  answered  are: 

•  Can  better  performance  be  obtained  using  the  original  compensator  but 
requiring  a  better  output? 

•  Can  the  compensator  be  reduced  and  the  resulting  output  be  as  good  or  better 
than  the  original  compensators? 

•  Can  a  better  output  be  obtained  by  increasing  the  order  of  the  compensator? 

When  the  Function  Minimization  subroutine  minimizes  the  function,  it  moves 

the  poles  and  zeroes  of  the  compensator  to  the  best  locations  in  the  S-plane.  If  the 
compensator  has  too  many  poles  and  zeroes,  the  program  tends  to  set  Z  =  P  for  the 
unneeded  poles  and  zeroes,  thus  the  compensator  order  may  be  reduced. 

In  many  cases,  even  though  the  poles  and  zeroes  may  not  exactly  cancel,  they  are 
placed  so  close  to  each  other  that  they  contribute  very  little  to  the  result  and  so  we 
may  efficiently  cancel  them.  Thus  one  can  start  the  design  by  simply  choosing  each  Kc 
to  have  numerous  zeroes  and  poles  as  can  be  seen  in  Table  7  and  Figure  3.2,  which  is 
the  original  pre-compensator.  However,  this  thesis  will  show  us  that  we  don’t  need 
that  many  poles  and  zeroes  for  the  choice  of  compensators. 

In  this  thesis,  as  shown  in  Table  9  and  Figure  3.4,  the  reduced  order  pre¬ 
compensator  was  used. 

The  optimum  parameter  values  of  this  reduced  order  pre-compensator,  which 
reduces  the  corresponding  cross-coupling,  can  be  determined  by  Function 
Minimization  and  Table  10  shows  the  optimum  parameter  values  of  the  reduced  order 
pre-compensator  in  Table  9  by  Function  Minimization. 

Figure  4.2  show's  the  positive  damping  effect  on  the  system’s  time  response  when 
optimum  parameter  values  are  used  in  the  reduced  pre-compensator. 

Figures  4.3  through  Figure  4.5  show  the  difference  between  Roberts'  method  and 
the  Function  Minimization  method,  where  the  continuous  line  is  the  output  of  Roberts 
and  the  dotted  line  is  the  output  of  Function  Minimization. 

As  shown  in  Figure  4.4  and  Figure  4.5,  cross-coupling  is  decreased  dramatically 
in  the  case  of  Function  Minimization. 

C.  CHOICE  OF  DESIRED  OUTPUTS 

One  can  always  choose  the  ideal  output  as  the  desired  output  for  the  Function 
Minimization  subroutine.  However,  this  may  not  be  a  good  choice  because  that 
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TABLE  10 

OPTIMUM  PARAMETER  VALUES  OF  REDUCED  PRE¬ 
COMPENSATOR 


RKn( S) 

1 

8.  6S  +  1 

RK12(S) 

K.2(-72.14S3  +  1.98925S2  -  2.  3S 

76.  2SJ  +  28.  63  S2  +  11.  5S  +  1 

RK^(  S) 

= 

1 

rk31(S) 

. 63525 

32. Q5S2  +  2.  975S  +  1 

RK32(  s  ) 

K.3(  1741.  25S:  -  1.  94375S  +  1) 

9.25S3  +  111.25S2  +  47.  IS  +  1 

rk33(S) 

1 

11.  4S  +  1 

RK13(S) 

2 

RK^S)  =  RK33(  S)  =  0 

particular  output  may  be  impossible  to  achieve.  Although  the  minimization  process 
will  determine  a  closest  fit,  if  the  cost  function  is  a  "least  squares"  function  the 
solution  may  not  be  acceptable.  For  example,  with  the  ship  problem  of  this  thesis,  the 
ideal  output  might  be 

•  Turn  radius  of  two  ship  lengths 

•  Zero  speed  change 

•  Zero  roll  angle 

None  of  these  characteristics  are  possible.  The  results  of  a  Function 
Minimization  design  will  not  satisfy  any  of  them  and  the  design  achieved  may  not  be 
acceptable.  The  desired  output  must  therefore  be  chosen  realistically,  i.e.,  -within  the 
physical  capabilities  of  the  system.  For  the  ship  control  problem,  outputs  were  chosen 
based  upon  the  results  obtained  by  Roberts.  In  order  to  obtain  better  performance  the 
desired  outputs  were  chosen  to  be  similar. 
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Figure  4.3  Comparison  between  Roberts  and  F.M  at  Yaw  output. 
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-  ISO 


Figure  4.4 
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Comparison  between  Roberts  and  F.M  at  Roll  output. 
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Figure  4.5  Comparison  between  Roberts  and  F.M  at  Speed  output. 
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This  thesis  chose  3  desired  system  outputs  ,  each  consisting  of  a  specification  for 
heading,  roll  and  speed.  These  desired  output  curves  (AFGEN1,  2  and  3)  were 
obtained  by  modifying  the  known  outputs  (Figure  3.3)  using  Arbitrary  Function 
Generators  (AFGEN)  provided  by  DSL. 

Figure  4.6,  Figure  4.7  and  Figure  4.8  show  AFGEN1,  AFGEN2  and  AFGEN3 
respectively. 

We  can  know  that  the  AFGEN 3  is  better  than  the  other  two  cases  by  the  trial 
and  error  method.  Therefore  this  thesis  will  use  the  AFGEN3  for  desired  outputs  for 
another  simulation. 

D.  GOST  FUNCTIONS 

When  using  Function  Minimization  as  a  design  tool,  a  cost  function  must  be 
chosen.  This  cost  junction  is  usually  3n  integral. 

One  oossibie  procedure  s  to  ..noose  'he  desired  performance  as  a  reference  ina 
>elect  i  cost  -'unction  wmen  is  "he  integral  of  'he  square  of  'he  difference  between 
uesireu  output  and  actual  output.  When  the  system  has  several  outputs,  the  cost 
function  must  consider  ail  of  mem.  usually  as  a  weighted  sum. 

E.  WEIGHTING  FACTORS 

For  the  mree  outputs  system,  the  cost  function  is  of  the  form 
E  -  X , J E J 2  dt  «-  X,jE;2  dt  -  X3jE32  dt. 

There  are  no  fixed  rules  for  choosing  the  values  factors  X,,  X,  and  Xr  One 
approach  is  as  follows: 

•  Select  the  output  whicn  is  considered  most  important,  say  n  1.  Use  this  as  a 
reference  and  -.elect  X,  =  1.0. 

•  Base  the  vaiue  assigned  to  the  second  weighting  factor.  X,  in  'his  case,  in  'he 
.importance  a:  output  =  I  wtt.n  respect  to  output  a  i.  If  output  s  2  is  equaily 
important  as  output  a  1.  'hen  X,  =  1.0.  If  it  is  naif  as  important,  then  X,  = 


•  The  third  weighting  factor  is  chosen  in  like  manner.  If  the  third  output  is  less 
important  than  the  first,  then  X.  <  X.  If  it  is  also  less  important  than  the 
second  output,  'hen  X,  X. 

The  actual  numbers  chosen,  for  example  X,  =  0.1.  are  simply  estimates  based  on 
experience.  The  designer  may  decide  to  change  'hem  after  studying  simulation  results 
of  a  first  design. 

In  the  case  of  ship  control,  the  primary  output  is  heading.  Therefore,  assign  Xh  a 
value  of  1.0.  Some  change  in  speed  is  unavoidable  and  the  desired  response  should 


« 
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show  this.  However,  if  the  desired  response  cannot  be  achieved  perhaps  more 
fluctuation  in  speed  would  be  acceptable.  For  example,  a  value  of  approximately  0.1 
could  be  assigned  to  If  the  desired  roll  output  is  twice  as  important  as  the  speed 
output,  then  would  be  assigned  a  value  of  0.2. 

Table  11  shows  each  case  for  simulation  using  weighting  factors. 


TABLE  11 

SIMULATION  CASES  USING  WEIGHTING  FACTORS! W.F.) 


CASE  :  A_  !  A  j  X  !  OUTPUTS 

1  n  r  I  s 


f  ! 

.2  .1  Figure  4.  9 

.6  .1  Figure  4.  10 


CASE4  is  from  Table  10,  which  uses  Xh  =  1,  X,f  =  1  and  -  1  and  this  was 
repeated  in  Table  1 1  for  comparison  when  we  consider  the  weighting  factors. 

Table  12  through  Table  15  show  the  elements  of  the  resulting  compensators  for 
each  of  the  above  when  we  used  the  F.M.  approach. 

When  F.M.  was  used  for  ail  coefficients  in  compensators.  RKp,  RK3l  and  RKp 
were  the  parameters  that  the  F.M.  subroutine  adjusted. 
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Figure  4.11  Step  Response  when  Xh  =  l, 
at  Rudder  Demand. 


,„**  .6  and  X  =  .4 

I  5 
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TABLE  12 

OPTIMUM  PARAMETER  VALUES  OF  REDUCED  PRE¬ 
COMPENSATOR  AT  CASE  1  W.E. 


RK„(S) 

RK12(S) 
RK22(  S) 
rk31(S) 

RK32(S) 

rk33(S) 

RKj3(S) 


1 

8.  6S  +  1 

Kr,(  -47.  64S'  +  2.  0C0625S2  -  2.  3S  +  1) 
67.  7S3  +  27._33S2  +  11.  375S  +  1 

1 

. 577125 

20.  425S2  +_  1.  8625S  +  1  ~ 

KcJ(  1915.  625S2  -  2.  0C1875S  +  1) 

"90.  625S3  +  53,  125S2  +  23.85S  +  1 

_ 1 

11.  4S  +  1  ” 

RKjjlS)  =  RK33(S)  =  0 


TABLF.  13 


OPTIMUM  PARAMETER  VALLES  OE  REDUCED  PRE¬ 
COMPENSATOR  AT  CASE  2  W.F. 


i 
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TABLE  14 

OPTIMUM  PARAMETER  VALUES  OF  REDUCED  PRE¬ 
COMPENSATOR  AT  CASE  3  W.F. 


RKn(S)  = 


8.  6S  +  1 


K,(-56.14S3  +  1.996S2  -  2.  3S  +  1) 


TABLE  15 

OPTIMUM  PARAMETER  VALUES  OF  REDUCED  PRE¬ 
COMPENSATOR  AT  CASE  4  W.F. 


RKU(S) 

RK12(S) 


8.  6S  +  1 

K,,(-72.14S3  +  1.  98925S2  -  2.  3S  *  1) 
76.  2S3  +  28.  63  S2  ♦  11.  5S  +  1 


R^S)  =  1 


RK31(S) 

rk32(s) 


. 63525 

32. 05S2  +  2.  975S  +  1  ' 

Kc.(  1741.  25S2  -  1.94375S  +  1) 
9.25S3  +  111.25S2  +  47.  IS  +  1 
1 


RK33(S)  = 


11.  4S  +  1 


RK13(S)  =  RK21(S)  =  RK33(S)  =  0 


From  simulation,  it  appears  that  the  most  desirable  output  can  be  obtained  at  - 
Xh  =  I  ,  Xf  =  .2  and  Xs  =  .1. 

Figure  4.13  through  Figure  4.15  show  the  comparison  between  the  Figure  3.4  for 
Roberts,  Figure  4.2  for  CASE  4  and  Figure  4.9  for  CASE  L  Where  the  output  of 
Roberts  is  from  Table  9,  which  gives  the  final  reduced  pre-compensator  elements,  the 
output  of  CASE  4  is  from  Table  10,  which  gives  the  optimum  parameter  values  of 
reduced  order  pre-compensator  by  Function  Minmization.  The  output  of  CASE  1  is 
the  best  case  when  we  consider  the  weighting  parameters  for  CASE  1,  2  and  3  in  Table 
11. 

The  continuous  line  is  used  for  Roberts,  the  dotted  line  for  CASE  4  and  the 
dashed  line  for  CASE  1  outputs. 

Figure  4.16  through  Figure  4.18  show  the  comparison  between  the  above  three 
cases  for  Fin,  Rudder  and  Power  respectively. 
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Figure  4.15  Comparison  between  I  ig  3.4,  Fig  4.2  anti 

at  Speed  Output. 
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Figure  4.16  Comparison  between  Roberts,  CASE  4  and  CASE  I 
for  Fin  Output  at  3°  Rudder  Demand. 
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Figure  4.17  Comparison  between  Roberts,  CASE  4  and  CASE  1 
for  Rudder  Output  at  3°  Rudder  Demand. 
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Figure  4.18  Comparison  between  Roberts,  CASE  4  and  CASE  1 
for  Power  Output  at  3°  Rudder  Demand. 


The  simulation  obtained  so  far  was  performed  using  three  degrees  for  rudder 
demand  and  we  can  know  that  there  is  no  problem  about  maximum  values  of  input 
and  output  shown  by  Table  3  for  each  output. 

As  can  be  seen  in  Figure  4.16  through  Figure  4.18,  CASE  1  is  better  than  the 
other  two  cases.  But  we  have  to  consider  carefully  the  physical  realization  in  Figure 
4. IS  because  the  power  curve  varies  very  fast.  Therefore  when  we  design  the  pre¬ 
compensator  using  weighting  factors  or  other  factors,  we  have  to  always  consider 
physical  realization  limitations. 

Figure  4.19  through  Figure  4.21  show  the  comparison  between  the  above  three 
cases  for  Fin,  Rudder  and  Power  when  we  use  twenty  degrees  for  rudder  demand. 

From  the  above  outputs,  we  can  see  that  CASE  4  is  better  than  the  other  two 
cases  when  we  use  twenty  degrees  for  rudder  demand. 

F.  VARIATION  OF  NUMBER  OF  POLES  AND  ZEROES 

As  shown  by  Table  7,  the  original  pre-compensator  by  Roberts  requires  complex 
mathematical  equations  for  a  complete  and  detailed  description.  For  many  problems 
in  the  analysis  and  design  of  many  modem  dynamic  systems,  a  simplified  description, 
i.e.,  a  low  order  model,  is  adequate  and  desirable.  This  thesis  is  concerned  with  the 
development  of  such  low  order  models  for  the  pre-compensator. 

Two  types  of  situations  are  commonly  encountered  in  practice: 

•  A  system  exists  and  can  be  tested,  but  its  equations  are  not  well  known  or  not 
clearly  defmed. 

•  A  high  order  complex  model  of  the  system  is  known  and  can  be  used,  but 

it  is  undesirable  for  the  problem  to  be  studied. 

In  either  case  the  response  of  the  system  to  a  chosen  signal  can  be  obtained,  and 
a  low  order  model  developed  which  has  essentially  identical  outputs  for  identical  inputs 
as  the  higher  model. 

By  carefully  planned  studies  it  should  be  possible  to  see  how  much  reduction  in 
order  can  be  achieved,  how  closely  the  behavior  of  the  low  order  compares  to  that  of 
the  system  and  perhaps  a  best  or  optimum  order  can  be  found  for  the  reduced  order 
models. 

In  addition,  by  careful  selection  and  classification  of  the  pole-zero  geometry  of 
the  high  order  system,  it  is  hoped  that  a  correlation  may  be  found  between  such 
geometry  and  that  of  the  best  low  order  model.  In  any  event,  it  is  anticipated  that 
some  rules  may  be  established  for  the  choice  of  the  number  of  poles  and  zeroes  in  the 
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Figure  4.21  Comparison  between  Rohcrts,  CASE  4  and  CASE  l 
for  Power  Output  at  20°  Rudder  Demand. 


51 


•  YslOX'A*  jl£lYO. 


low  order  model,  such  that  evaluation  of  the  pole,  zero  and  gain  will  require  only  a  few 
computer  runs. 

As  shown  in  Figure  A.l  of  System  Block  Diagram  for  Simulation,  Kn  and  Kp 
afFect  the  Fin  output  of  the  system  and  K31,  K32  and  K33  affect  the  Power  output  of 
the  system.  Therefore  when  we  change  the  order  of  these  elements,  outputs,  i.e. ,  roll 
and  speed  will  be  changed  by  them  also. 

Unfortunately  there  is  no  known  mathematical  basis  for  choosing  the  "best" 
order  for  a  low  order  model  for  each  element  of  the  pre-compensator.  For 
convenience,  this  thesis  has  changed  the  order  of  Kp  and  K33  for  two  cases  only. 

Table  16  and  Table  17  show  the  values  of  element  for  two  cases  and  Figure  72 
through  Figure  25  show  the  outputs  for  comparison  between  original  reduced  order 
pre-compensator  by  Table  10  and  the  above  two  cases  respectively.  Other  elements  are 
the  same  as  in  Table  10. 


TABLE  16 

OPTIMUM  PARAMETER  VALUES  OF  REDUCED  PRE- 

COMPENSATOR 

WHEN  WE  CHANGE  ONLY 

KU(S) 

97. 18759S  +  3. 94875 

7.53125S2  +  21.075S  +  12.  025 

K12(S) 

Kc2(-72.14S3  +  1.  98925S2  -  2.  3S  +  1) 

76.  2S3  +  28.  63S2  +  11.  5S  +  1 

K:,(  S  ) 

= 

1 

K3l(S) 

. 63525 

32. 05S2  +  2.  975S  +  1 

k32(S) 

Kc3(  1741.  25S2  -  1.  94375S  +  1) 

9.  25S3  +  111.25S2  +  47.  IS  +  1 

k33(S) 

1 

11.  4S  +  1 

K13(S) 

"SZ 

K21(S)  =  K33(S)  =  0 
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By  the  same  algorithm,  we  can  approach  the  appropriate  order  of  pre¬ 
compensator  and  the  values  of  each  element  for  better  design. 

As  can  be  seen  in  Figure  4.22  through  Figure  4.25,  when  we  put  the  additional 
zero  and  pole  to  the  element  of  the  pre-compensator,  the  cross-coupling  of  the  roil 
output  is  decreased  dramatically. 
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Figure  4.22  Comparison  between  Table  10  and  Table  16 
for  Roll  Output. 
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Figure  4.23  Comparison  between  Table  10  and  Table  16 
for  Speed  Output. 
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Figure  4.25  Comparison  between  Table  10  and  Table  17 
for  Speed  Output. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

In  this  thesis,  the  development  of  a  pre-compensator  to  reduce  the  undesired  and 
highly  noniinear  cross-coupiing  effects  in  the  maneuvering  characteristics  of  a  modern 
warship  has  been  presented.  Simulation  results  have  shown  that  Function 
Minimization  procedures  for  coordinated  steering  of  a  surface  '•nip  would  significantly 
improve  ship  stabiiitv.  minimize  loss  of  soeed  and  reduce  the  interaction  intensity  of 
the  compensated  system. 

Simulation  resuits  have  also  shown  that  the  reduced  pre-compensator  can  be 
determined  by  Function  Minimization  directiy  from  an  ideai  pre-compensator  based  on 
given  'pecuicauons. 

B.  RECOMMENDATIONS 

Computer  simulation  for  Function  Minimization  leads  to  the  following 
recommendations: 

•  This  thesis  uses  the  "HOOKE"  subroutine  of  DSL  for  Function  Minimization. 
It  has  certain  constraint  parameters,  i.e.,  UMAX,  CFT.  Simulation  outputs 
are  changed  by  changing  initial  values  and  the  step  size  for  Function 
Minimization.  Therefore,  for  given  specifications  initial  values  and  step  size  for 
Function  Minimization  by  constraint  parameters  should  be  determined  by 
experience. 

•  In  rhis  thesis,  it  has  been  shown  that  Function  Minimization  can  be  used  for 
various  cases  as  shown  in  Table  11.  A  particular  case  for  needed  specification 
should  be  determined  and  simulated  by  trial  and  error  method  for  element 
optimum  values  which  have  undesired  cross-coupiing  effects. 

•  In  this  thesis,  it  has  been  assumed  that  the  stabilizer  fm  to  yaw  cross-coupling 
term,  G^fS),  is  a  null  entry'  for  the  class  of  warship  considered.  However,  it  is 
recommended  that  cross-coupling  between  these  parameters  be  considered  in 
further  studies  on  general  surface  ships. 

•  Further  research  should  investigate  ship  characteristic  constraints. 

•  Further  research  should  investigate  the  effects  of  various  sea  state  conditions, 
speeds  and  maneuvering. 
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APPENDIX  A 

SYSTEM  BLOCK  DIAGRAM  FOR  SIMULATION 


Figure  A.l  System  Block  Diagram  for  Simulation. 
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APPENDIX  B 

CONSTRAINT  PARAMETERS  OF  FUNCTION  MINIMIZATION 

To  use  HOOKE,  the  user  must  initialize  the  following  arguments  and,  in  the 
main  program. 

CALL  HOOKE  (X,  STEP,  N,  ITMAX,  CFTOL,  ALPHA,  BETA,  CF,  Q,  QQ, 
W,  IPRINT,  MINMAX) 

Ail  of  these  arguments  must  be  initialized  in  MAIN,  except  for  X,  CF,  Q,  QQ, 
and  W.  Recommended  values  for  ALPHA  and  BETA  are  ALPHA  =  2..  BETA  =  0.5. 
All  of  the  arrays,  i.e.,  X,  STEP,  Q,  QQ  and  W,  must  be  declared  and  dimensioned  in 


TABLE  18 

THE  LIST  OF  PARAMETERS 

ARGUMENT 

MEANING 

X 

the  array  of  N  parameter  values.  The  user  must 
suoply  the  initial  guesses,  either  in  the  DSL 
program  or  in  MAIN. 

STEP 

an  array  of  dimension  N  containing  the  initial 
stepsizes  to  be  used  in  the  seared. 

N 

the  number  of  oarameters  ( a  positive  integer, 
at  most  15). 

the  maximum  number  of  function  calls  no  be 
ITMAX  performed. 


the  error  in  the  criterion  function  to  be 
CFTOL  reached  before  the  orogram  terminates 

,( difference  between' the  current  value  and  the 
previous  stage  value). 


the  factor  of  (Y  -  X)  which  is  added  to  Y  to  cret 
ALPHA  IXNEW;  a  number  greater  than  or  equal  to  1. 

! the  stepsize  reduction  factor;  a  number  between 
BETA  !0.  and  1. 


CF  the  /alue  of  the  criterion  function. 


i 


Q/ QQ, w 


arrays  of  dimension  N,  to  be  used  as  work  SDace. 
Thev  must  be  declared  and  dimensioned  in  the 
MAIN  orogram. 


: an  integer  fiag:=  0  for  no  intermediate  printout 
IPRINT  |  ~  1  for  intermediate  printout  I 

I  of  X,  CF,  the  number  of  function  evaluation  and 
l notification  of  step- reduction. 


an  integer  flag:  =  -1  searches  for  a  minimum 

=  +1  searches  for  a  maximum 


MINMAX 


APPENDIX  C 

COMPUTER  PROGRAM  FOR  UNCOMPENSATED  SYSTEM 

TITLE  SIMULATION  OF  UNCOMPENSATED  SYSTEM 
ARRAY  Al(  1),B1(3),A2(2),B2(4),A3(  1),B3(5),A4(2),B4(4),... 

A5(  1  ),B5(3).A6(1  ),B6(2) 

TABLE  AKII-  l.Bl(  l-3)»4..24.1.A2(l-2)  = -8.57, l.B2(l-4)»  53.3.17. r.... 

9.52.1,A3{’1)  =  l.B3(l-5)  =  12.32.25.1 1.2.1.0.A4(I-2)  =  10.1 .... 

B4(l-4)  =  240,5S,26.1.A5(1)=  1.35(1-3)=  240.34.  i.Ab(  1)=  1.... 

B6<  1-7)=  24.1 

CONST  K  !  1  =  0. 1 14.K  12  =  4. 1 8.K22  =  .01  ,K3 1  =  0.058.  K32  =  0.096.  K3  3  =  u.  1,... 

RL  =  0.,YW  =  0..SP  =  0. 

DERIVATIVE 

FIN  =  0.0*STEP(0) 

RUDDER®  3.'!STEP(0) 

POWER  =  0/'STEP(0) 

ROLLl  =  TRNFR(0.2.RL.AI.B1.K  1 1  *  F  T  N ) 
ROLL2-TRNFR(1.3.RL,A2.B2.Ki:*RUDDER) 

ROLL=  ROLL!  -  ROLL2 

YAW  =  TRNFR(0,4,YW.A3,B3.K22*RUDDER) 

SPEED  1  =  TRNFR(  1,3.SP,A4,B4,K31*FIN) 

SPEED2  =  TRNFR(0,2.SP.A5,B5,K  32*  RUDDER) 

SPEED3  =  TRNFR(0,l,SP,A6,B6.K33*POWER) 

SPEED®  -SPEED1-SPEED2  +  SPEED3 
CONTROL  FINTIM  =  80. 

*R1NT  L. ROLL. YAW, SPEED 
SAVE  (SI)  0.1. ROLL, YAW, SPEED 

GRAPHfGl  S1.DE=TEK618,PO  =  0,.5)  TIMEILE  =  S..UN=  SEC)  ROLL(LO  =  -l.LI  =  1.... 
SC  =  .25.NI  =  S),YAW(  LO  =  -l.LI  =  3. SC  =  .25, NT  =  5).... 

SPEED!  LO=  -l.LI  =  4.SC  =  25.N1  =  S) 

LABEL).  GUSTEP  RESPONSE  FOR  UNCOMPENSATED  WARSHIP  MODEL  AT  12  KTS 

LABEL) Gl  TO  %  RUDDER  DEMAND 

END 

STOP 
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APPENDIX  D 

COMPUTER  PROGRAM  FOR  ORIGINAL  PRE-COMPENSATOR 

TITLE  SIMULATION  OF  ORIGINAL  PRE-COMPENSATOR 
ARRAY  A1(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4(4),... 
A5(1),B5(3),A6(1),B6(2) 

ARRAY  Cl(  1),D1(4),C2(4),D2(4),C3(2),D3(6).C4(6),D4(7),C5(1),D5(4) 

TABLE  Al(l)=  l,Bl(l-3)  =  4,.24,l,A2(l-2)  =  -8.57,l,B2(l-4)  =  53.3,17.17,... 
9.52,1,A3(  1)  =  1,B3(  1-5)=  12,32.25,1 1.2, 1,0,  A4(  1-2)  =  10.1,... 

B4(  1  -4)  =  240,58,26, 1  ,A5(  1 )  =  l,B5(l-3)=  240, 34,1,  A6(l)=  1,... 

B6(  1-2)=  24.1 

TABLE  Cl(l)=  1,D  1(1-4)=  53.19,17.136,9.5,1,... 

C2(l-4)=  -34.143,2,-8.29,1,... 

D2(l-4)  =  53.19,17. 13, 9.5,1, C3(  1-2)=  10,1,... 

D3(  I-6)=  534,278.182.6.46.15,1 1.5,1.... 

C41 1  -6)  =  80050. 15650,20665.2255,3, 1 .... 

D4(  1-7)  =  5352, 3326,21 13,646, 162,21. 6,1,... 

C5(l)=  l,D5(T-4)=  100,30,12,1 

CONST  K1 1  =  0. 1 14, K 12  =  0. 18, K22  =  .0l, K3I  =  0.058, K32  =  0.096, K33  =  0.1,... 

RL  =  0.,YW  =  0.,SP  =  0. 

CONST  GC 10  =  0,GC20  =  0,GC30  =  0.GC40  =  0,GC50  =  0 
DERIVATIVE 

RIN  =  0.0*STEP(0) 

YIN  =  3.*STEP(0) 

SIN  =  0.*STEP(0) 

El  =  RIN-ROLL 
E2  =  YIN-YAW 
E3  =  SIN-SPEED 

FIN1  =  TRNFR(0,3,GC10,C1,D1,E1) 

FIN2  =  TRNFR(3,3,GC20,C2,D2,-I.6*E2) 

FIN  =  FIN1  +  FIN2 

ROLL1  =  TRNFR(0,2.RL,A1,B1,K11*FIN) 

RUDDER®  YIN-YAW 

ROLL2  =  TRNFR(  1.3, RL,A2,B2,K  12*  RUDDER) 

ROLL®  ROLL1  +  ROLL2 

YAW  =  TRNFR(0,4,YW.A3,B3,K22*RUDDER) 

SP1  =  TRNFR(  1,5,GC30,C3,D3,0.579:!'E1) 

SP2  =  TRNFR(5,6,GC40,C4,D4,0.0456*E2) 

SP3  =  TRNFR(0,3,GC50,C5,D5,E3)  • 

SPEED1  =  TRNFR(1,3,SP,A4,B4,K31*FIN) 

SPEED2  =  TRN  F  R(0, 2, SP,A5,B5,K32*  RUDDER) 

POWER®  SP3-SP2-SP1 

SPEED3  =  TRNFR(0,1,SP,A6,B6,K33*POWER) 

SPEED®  SPEED3-SPEED2-SPEED1 
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CONTROL  FINTIM  =  80. 

•RINT  1., ROLL, YAW, SPEED 
SAVE  (SI)  0.1, ROLL, YAW, SPEED 

GRAPH(G1/S1,DE  =  TEK61 8,PO  =  0,.5)  TIME(LE  =  8., UN- SEC)  ROLL(LO  =  -2, LI  =  1,... 
SC  =  .5, NI  =  8),YAW(  LO  =  -2, LI  =  3, SC  =  .5, NT  =  8),... 

SPEED(LO  =  -2, LI  =  4, SC  =  .5, NT  =  8) 

LABEL(G1)STEP  RESPONSE  FOR  COMPENSATED  WARSHIP  MODEL  AT  12  KTS 

LABEL(G1)10  %  RUDDER  DEMAND 

END 

STOP 


APPENDIX  E 

COMPUTER  PROGRAM  FOR  ORIGINAL  REDUCED  PRE¬ 
COMPENSATOR 

TITLE  SIMULATION  FOR  ORIGINAL  REDUCED  ORDER  PRE-COMPENSATOR 
ARRAY  A1(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4I4),... 

A5(I).B5(3).A6(I  ),B6(2) 

ARRAY  C  1(  1  ),D 1('2).C2(4).D2(J),C3(  1  ).D3( 3LC-I3 ),D4( 4'i,C5('  1  i.D5<  2 1 
TABLE  AH1)=  l.BHl-3)- J..24.1.A2(I-2)  =  -3.57.I.32(  l-4)  =  53.3.17.1".... 

9.52.i.A3(  1  )=  l.B3n-5)=  12.32.25.1 1.2.1.0.A4U-2)  =  10.1.... 

34(1-4)=  240.5S.26,i.A5(l)=  1.35(  1-3)=  240.34. 1.A61  1  )=  i,... 

B6(T-2)=  24.1 

TABLE  Cl(  1)=  l.D  1(  1-2)=  8.6.1.C21  1-4)=  -34. 14, 2.-8. 3.1.... 

D2(  1-4)=  53.2.1".  13. 9.5,1.  C3(  1)  =  I.D3d-3)  =  14.3.1.3.1.... 

C4(  1  -3 )  =  2000.-2. 1  .D4/1  -4)  =  ;  30.25. 1 2.6.  i .... 

C5(T‘=  1.D5G-2)  =  11.4.1 

CONST  K 1 1  =  0. 1 14. K 12  =  0. 1S.K22  =  .01. K2 1  =  0.05S.K32  =  0.096.K33  =  0.1.... 

RL  =  0..YW  =  0..SP  =  0. 

CONST  KC 1  =  -  1.6.KC2  =  0.579. KC3  =  0.0456,... 

GC 10  =  0.GC20  =  0.GC30  =  0.GC40  =  0.GC50  =  0 
DERIVATIVE 

RIN  =  0.0*STEP(0) 

YIN  =  3.*STEP(0) 

SIN  =  0.*STEP(0) 

El  =  RIN-ROLL 
E2=  YIN-YAW 
E3  =  SIN-SPEED 

PINT  =  TRNFR(0, 1.GCI0.C1.DI.E1) 

FIN2  =  TRNFR(3.3,GC20,C2,D2,KC1*E2) 

FIN  =  FIN  1  -  FIN2 

ROLL!  =  TRNFR(0.2.RL.AI.31.K1 1*FIN) 

RUDDER  =  YIN-YAW 

R0LL2=TRNFR(1,3.RL.A2,B2,K12*RUDDER) 

ROLL=  ROLL1  +  ROLL2 

YAW=  TRNFR(0.4.YW.A3.B3.K  22*  RUDDER) 

SPl  =  TRNFR(0.2.GC30.C3.D3.KC2*EI) 

SP2  =  TRNFR(2.3,GC40.C4.D4,KC3a‘E2) 

SP3  =  TRNFR(0,1,GC50,C5,D5.E3) 

SPEED1  =  TRNFR(  1,3.SP,A4,B4,K31SFIN) 

SPEED2  =  TRNFR(0,2,SP,A5,B5,K32*RUDDER) 

POWER  =  SP3-SP2-SP1 

SPEED3  =  TRNFR(0,1,SP,A6,B6,K33*POWER) 

SPEED  =  SPF.ED3-S  PEED  2-S  PEE  D1 
CONTROL  FINTIM  =  80. 
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*RINT  1., ROLL, YAW, SPEED 
*RINT  1., FIN, RUDDER, POWER 
SAVE  (SI)  0.1, ROLL.YAW, SPEED 

GRAPH(G  1/S  1  ,DE  =  TEK6 1 8,PO  =  0,.5)  TIME(LE  =  8., UN- SEC)  ROLL(LO=  -2, LI  =  1,... 
SC  =  .5,NI  =  8),YAW(  LO  =  -2, LI  =  3.SC  =  .5,NI  =  8),... 

SPEED(LO  =  -2,LI  =  4, SC  =  .5,NI  =  8) 

LABEL(G1)STEP  RESPONSE  FOR  COMPENSATED  WARSHIP  MODEL  AT  12  KTS 

LABEL(G1)FIGL'R£  3.4 

END 

STOP 


APPENDIX  F 

COMPUTER  PROGRAM  FOR  AFGEN  SUBROUTINE 

TITLE  SIMULATION  FOR  desired  outputs  using  afgen  subroutine 
AFGEN  RL  =  0,0,1, -.03, 4,-.  1,10,0,16,0,23,0,29,0,36,0,... 
42,0,48,0,54.0,61,0,67,0,80,0 

AFGEN  YW  =  0,0, 1 ,. 00006, 4, . 006, 8, . 04, 10,. 07, 20..29, 30, .55, 40, .8, 80, 1.61 
AFGEN  SD  =  0,0, 1  ,-.004. 10, -.0 1 ,2 1 ,0,3 1 ,0,4 1 ,0,5 1 ,0,70,0,80,0 
DERIVATIVE 

X  =  TIME 

ROLL  =  NLFGEN7RL.X) 

YAW  =  NLFGENYYW.X) 

SPEED  =  NLFGEX(SD,X) 

CONTROL  FIN'TIM  =  80 
*RINT  1. ROLL.  YAW, SPEED 
SAVE  iSl)  0. 1. ROLL. YAW. SPEED 

GRAPH1GI  S1,DE  =  TEK618.PO  =  0,.5)  TIME(LE=  8., UN  =  SEC), ROLL(LO 
SC  =  .5, XI  =  8).YAW<TO=  -2, LI  =  3,SC=  .5, NT  =  8),... 

SPEED(  LO  =  -2.LI  =  4, SC  =  .5, XI  =  8) 

LABEL(Gl)  AFGX  1  FOR  FIGURE  4.6 

END 

STOP 


APPENDIX  G 

COMPUTER  PROGRAM  FOR  FUNCTION  MINIMIZATION 
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D  COMMON’HANDJ/FLAG,ER,K  1,K2,K3,K4,K5,K6,K7,K8,K9,K  10,K1 1,K12 
D  COMMON  HAN'DJ  K13.K  14 

TITLE  SIMULATION  FOR  FUNCTION  MINIMIZATION  SUBROUTINE 
ARRAY  A1(1),B1(3),A2(2),B2(4),A3(1),B3(5),A4(2),B4(4),A5(1),B5(3),... 

A6(  1  ),B6<  2) 

STORAG  Cl(  l).Dl(2),C2(4),D2(4),C3(  1),D3(3),C4(3),D4(4),C5(1),D5(2) 

TABLE  A 1  ( I )  =  l.BK 1-3)  =  4,.24,1,A2(  1-2)=  -8.57,1,B2(  I -4)  =  53.3,17. 17.... 

9.52.1,A3(1)  =  I.B3(  1-51=  12,32.25.1 1.2, 1.0,  A4<  1-2)=  10,1,... 

B4(  1-4)  =240, 58,26,1, A5(l)=  l,B5(l-3)  =  240,34,1,... 

A6(  1 )  =  1,B6(  1-2)  =24,1 

CONST  KOI  =  -34.14.K20  =  2,K30  =  -8.3,K40=  53.2.K50=  17.13. K60  =  9.5, K70=  .579.... 

KS0=  14,8. K90=  1.3,K100  =  2000.KI  10= -2,K120=  130,K130=  25.... 

K 140  =  12.6 

CONST  01=1 14.C12  =  .  18,C22=  01,C31  =  .058,C32  =  .096.C33  =  .  1,IC  =  0,... 

KC1  =  -1.6.KC3  =  0-0456 

PA  RAM  K 1MIN  =  -50.KI MAX=  200, K2 MIN  =  .01.K2MAX  =  200.K3MIN  =  -10.K3MAX  =  200.... 
K4M  IN  =  .0l.K4.MAX  =  200,K5.MIN  =  .01,K5.MAX  =  200.K6MIN  =  .01.K6MAX  =  200... 
K7MIN  =  .01.K7MAX  =  200.K8M1N  =  .01,K8MAX=  200,K9MIN  =  .01,K9MAX  =  200,  .. 

K  I0.M  I N  =  .0  UK  10MAX  =  3000, K 1 1  MIN  -  - 10, K 1 1  MAX  =  200.K  1 2M I N  =  .0 1 ... . 

K 1 2. MAX  =  200.K13MIN  =  .01,K  l3.VIAX  =  200,K14MIN  =  .01,K  14MAX=  200 
AFGEN  RL  =  0.0,1, *  03, 4, -.1,1 0,0, 16. 0,23, 0,29, 0,36,0,... 

42,0,48,0,54,0,61,0,67,0,80,0 

AFGEN  YW  =  0,0,1. .00006, 4, .006, 8, .04, 10..07.20, .29, 30, .55,40, .8, 80,1. 61 
AFGEN  SD  =  0,0, 1.-.004, 10, -.01, 21, 0,31, 0,41, 0,51, 0,70, 0,80,0 
INITIAL  SEGMENT 

IRFLAG.LT.O.)  Kl  =  KOI 
IFi  FLAG.LT.0.)  K2=  K20 
IF(  FLAG.LT.0.)  K3=  K30 
IFIFLAG.LT.O.)  K4=  K40 
IF(FLAG.LT.O.)  K5=K50 
IF(FLAG.LT.O.)  K6=  K60 
IF(FLAG.LT.0.)  K‘=  K70 
IFIFLAG.LT. 0.)  K8=  K80 
I F( FLAG.LT.0.)  K9=  K90 
IF(FLAG.LT.O.)  K10=  K100 
I F( FLAG.LT.0.)  Kll  =  K110 
IFf FLAG.LT.0.)  K12=KI20 
IF(FLAG.LT  0.)  K13=  K130 
IF(FLAG  LT.O.)  K14=K140 
FLAG  =  FLAG  +  1 

IF((K1.LE.K1MIN).0R.(K1.GE.K1MAX))  THEN 
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K1  =  KOI 
ENDIF 

IF((K2.LE.K2MIN).OR.(K2.GE.K2MAX))  THEN 
K2=  K20 
ENDIF 

IF((K3.LE.K3MIN).OR.(K3.GE.K3MAX))  THEN 
K3  =  K30 
ENDIF 

IF((K4.LE.K4M IN).OR.(K4.GE.K4MAX))  THEN 
K4  =  K40 
ENDIF 

IFf(K5.LE.K5MIN).OR.i'K5.GE.K5MAX))  THEN 
K5-  K50 
ENDIF 

IF((K6.LE.K6MIN).OR.(K6.GE.K6MAX))  THEN 
K6=  K60 
ENDIF 

;f(K'.le.k:min..or.;ktge.k'max))  then 

K7»  K70 
ENDIF 

IF((  KS.LE.K8MIN).OR.(  KS.GE.K8MAX))  THEN 
KS  =  K30 
ENDIF 

I FO  K9.LE.K9MIN)  OR.(  K9.GE.K9MAX))  THEN 
K9=  K90 
ENDIF 

IF((K  10.LE-K  lOMIN).OR.(KlO.GE.KlOMAX))  THEN 
K10=  KlOO 
ENDIF 

IFKKI  1.LE.K1  IMIN).0R4K1 1.GE.K1 1MAX))  THEN 
K11-K1IO 
ENDIF 

I  Fi  (  K 1 2. LE. K 1 2M INT.OR.i  K 1 2.GE. K 1 2M AXV)  THEN 

k;:=  ki:o 

ENDIF 

IF((K13.LE.K13MIN).OR.(  K  13.GE.K  I3MAX))  THEN 
K 13  =  K 1 30 
ENDIF 

IFi'(KI4.LE.KI4MIN).OR.;K14.GE.KI4MAXi)  THEN 
K14-  K 140 
ENDIF 
CI(I)=  1 
Dl(  1)=  S.6 
Dl(2)=l 
C2(  1)  -  K 1 
C2(2)  =  K2 
C2(3)  =  K3 
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C2(4)=  1 
D2(l)=K4 
D2(2)=  K5 
D2(3)=K6 
D2(4)=  1 
C3(l)=  K7 
D3(  1)=  XS 
D3(2)=K9 
D3(3)=  l 
C4(  1 )  =  K10 
C4(2)=  Kll 
040)  =  1 
D4n )  =  :<  i : 

D4,'2)  =  K13 
D4(3)=  K14 
D4(4)  =  1 
C5(  I)  =  1 
D5(l)  =  11.4 
D5f  2)  =  l 
DERIVATIVE 

RIN  =  0.0*STEP{0) 

YIN  =  3.*STEP(0) 

SIN  =  0.*STEP(0> 

El  =  RIN-ROLL 

E2-  YIN-YAW 

E3  =  SIN-SPEED 

FIN1  =  TRNFR(0,1,IC,CI,D1,EI) 

FIN2  =  TRNFR(3.3,IC,C2,D2.KC1*E2) 

FIN  =  FIN1  -1-  FIN2 

ROLL1  =  TRNFR(0, 2, IC,A1,B1, CHAFIN’) 
RUDDER  =  YIN- YAW 

roll:  =  trnfri  i.3.ic.a:,32,ci:*rudder) 
ROLL  =  ROLL1  *  ROLL2 
YAW  =  TRNFR(0.4.IC.A3.B3.C22*RUDDER) 
SPI  =  TRNFR(0,2,IC,C3.D3.E1) 

SP2  =  TRNFR(2,3,IC,C4,D4,KC3*E2) 

SP3  =  TRNFR(O.I,IC,C5,D5,E3) 

SPEED1  =  TRNFR(1.3.IC,A4,B4.C31*FIN') 
SPEED2  =  TRNFR(0.2.IC.A5.B5.C32*RUDDER) 
POWER  =  SP3-SP2-SP1 

SPEED3  -  TRNFR(0,1,IC.A6.B6,C33*POWER) 
SPEED  =  SPEED3-SPEED2-SPEED1 
X  =  TIME 

YNL1  =  NLFGEN(RL.X) 

YNL2  =  NLFGEN(YW,X) 

YNL3  =  NLFGEN(SD.X) 

EE1  =  ((ROLL-YNLl)**2) 
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Ell  *  INTGRL(IC,EE1) 

EE2  =  ((YAW-YNL2)**2) 

E22=*  INTGRL(IC,EE2) 

EE3- ((SPEED- YNL3)**2) 

E33=  INTGRL(IC,EE3) 

E  =  El  1  +  E22+  E33 
CONTROL  FIN'TIM  =  50.DELT  =  .01 
TERMINAL 
ER  =  E 

END 

STOP 

FORTRAN 

IMPLICIT  REAL*8( A-H.O-Z) 

DIMENSION  Xi  Mt.sTEPf  14j,Q(  14).QQ(  14),W»14) 

STEP!  1)  =  4. 

STEP(  2)=  .001 
STEPf  3 )  = 

STEP(4)  =  4. 

STEPL5)  =  2. 

STEP(6)  =  1. 

STEPD  =  0.005 
STEP(  8)  =  1. 

STEP(9)  =  0. 1 
STEP(10)-  15. 

STEP(  11)=  .005 
STEP(12)  =  7. 

STEP(  13)=  5. 

STEP(  14)  =  2. 

N  =  14 

IT.MAX  =  200 
CFTOL-. 000000 1 
ALPHA  =  2 
BETA  =  .5 
IPRINT-0 
M  INMAX®  -1 

CALL  HOOKE(X, STEP, N\ITMAX,CFTOL,ALPHA, BETA. 
*CF.Q,QQ,W,IPRINT.MINMAX) 

STOP 

END 


APPENDIX  H 

COMPUTER  PROGRAM  FOR  COMPARISON  BETWEEN  ROBERTS 

AND  F.M. 


TITLE  SIMULATION  FOR  COMPARISON  BETWEEN  ROBERTS  AND  F.M. 
*****  ROBERTS  OUTPUT  ***** 

ARRAY  A  I(  I  ),B  l(3),A2( 2).B2(4),A3(  1  ),B3(5),A4(2),B4(4),... 

A5(  1  ),B5(3),A6(  1  ),B6(2) 

ARRAY  C  1(  l),Dl(2).C2(4),D2(4),C3(  r),D3(3),C4(3),D4(4),C5n),D5(2) 

TABLE  Al(  1)  =  l.Bli  1-3)  =  4,.24.1.A2(  1-2)  =  -8.57,1, B2(  1-4)  =  53.3. r.l“.... 

4.52.1 . A3f  D=  I.B3H-5)*  12.32.25.1 1.2.1.0.A4I 1-2)=  10.1,... 

B4(  1-4)=  240.5S.26.i,A5(l)=  l.B5(l-3)=  240.34.1,  A6(l)=  1,... 

B6\  1-2 1  =  24,1 

TABLE  C  1(1 )  =  l.Dl(l-2)  =  8. 6. 1  .C2(  1-4)  =  -34.14.2, -8.3.1... . 

D2(l-4j=  53.2.1".  13. 9.5.1. C3(  1  >=  l.D3(l-3)=  14.S.  1 .3.  i .... 

C4(  1-3  .  =  200'».-2.1.D41 1-4)  =  1 30, 25. 12.6,1.... 

C5(l)=  l.D5(l-2)=  11.4,1 

CONST  K 1 1  =  0. 1 14.K  12  =  0. 18.K22  =  .01.K31  =  0.058.K32  =  0.096, K33  =  0. 1 .... 

RL  =  0..YW  =  0.,SP  =  0. 

CONST  KC1  =  -1.6.KC2  =  0.579,KC3  =  0.0456,... 

GC 10  =  0.GC20  =  0.GC30  =  0.GC40  =  0.GC50  =  0 
DERIVATIVE 

RIN  =  0.0*STEP(0) 

YIN  =  3.  *STEP(0) 

SIN  =  0.*STEP(0) 

El  =  RIN-ROLL 
E2=  YIN-YAW 
E3  =  SIN-SPEED 

FIN1  =  TRNFR(0.l.GC10,Cl.Dl,El) 

FIN2  =  TRNFRl  3.3.GC20.C2.D2.KCI  *E2) 

FIN  =  FIN  1  -  FIN2 

ROLL1  =  TRNFR(0.2.RL,A1.B1,KI1*FIN) 

RUDDER  =  YIN-YAW' 

ROLL2  =  TRNFR(  1.3,  RL,A2,B2,K  12*  RUDDER) 

ROLL=  ROLL1  -+-  ROLL2 

YAW  =  TRNFR(0,4,YW.A3.B3.K22*RUDDER) 

SP1  =  TRNFR(0.2.GC30,C3.D3.KC2*E1) 

SP2  =  TRNFR(2,3,GC40.C4.D4,KC3*E2) 

SP3  =  TRNFR(0,l,GC50,C5.D5,E3) 

SPEED1  =  TRNFR(1,3.SP,A4,B4,K31*FIN) 

SPEED2  =  TRNFR(0.2,SP,A5,B5,K32*RU  DDER) 

POWER=  SP3-SP2-SP1 

SPEED3  =  TRNFR(0,1,SP,A6,B6,K33*POW'ER) 

SPEED  =  SPEED3-SPEED2-SPEEDI 
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*****  FUNCTION  MINIMIZATION  OUTPUT  ***** 

ARRAY  AAl(l)tBBl(3)AA2(2),BB2(4),AA3(l),BB3(5),AA4(2),BB4(4)t... 

AA5(1),BB5(3),AA6(1).BB6(2) 

ARRAY  CC1(I),DD1(2),CC2(4),DD2(4),CC3(  1),DD3(3) 

ARRAY  CC4(3),DD4(4),CC5(  1),DD5(2) 

TABLE  AA  1(1)=  l.BBK  1-3)  =  4, .24.1, AA2(  1-2)  =  -8.57,1,... 

BB2I  1-4)  =  53.3.17. P.9. 52.1.... 

AA3(U  =  1,BB3(  1-5)=  12,32.25.1 1.2.1.0.AA4(l-2)  =  10.1.... 

BB4(  1-4)  =  240.58.26. 1  ,AA5(  1 )  =  1  ,BB5(  1  -3)  =  240.34. 1  ,AA6(  1 )  =  1 .... 

3B6<  1-2)=  24,1 

TABLE  <:c; 1  ,=  LOOK  1-2)  =  3.b.l.CC2(l-4)  =  -72.14.l. 98925.-2.3.1  .. 

OD2  I-4i-  “6.2.28.03. 11. 5.i. DD3(l-3)  =  32.05.2.975.1.... 

» l-'i  =  i*41.25.-i.94375.1.DD4U-4)-'),25.1 11.25.47. 1.1.... 

CT  1  .  =  ..OD5I  1-2)=  !  1.4.1.CC2i  1)=  o3525 
CONS  r  K K  1 1  =  >».  1 1  J.KK  12  =  o.  1S.KK22*  .0I.KK31  =  0.058.KK32  =  0.096.KK33 
RRL  =  0..YYW  =  0..SSP  =  0..KKC1  =  -1.6.KKC3  =  0.0456 
'  'ON  ST  GGC ! •< -  <  UK3C20  -  o.OGC'0  -  O.GGC  JO  =  O.GGC50  *  • > 


R 1 N  -=•*.>  STE Pf r*  > 
.  TIN  =  '.  'STEP")  . 


SSIN  = ').  "STEP'Oi 
EE  1  =  RR1N-RROLL 
EE2  =  YYIN-YY.AW 
EE3  =  SS1N-SSPEED 

FFIN1  =  TRNFRi  0. 1 .GGC10.CC l.DD LEE i  > 

FFIN2  =  TRNFRi  3.3.GGC20.CC2.DD2.KKC1  *EE2) 

FFIN-  FFIM-FFIN2 

RROLL1  =  TRNFRo).2.RRL.AAl.BBl.KKl  1*FF1N) 

RRL  DDE  =  YYIN-YYAW 

RROLL2  =  TRNFRi  1.3.RRL.AA2.BB2.KK  12*RRU'DDE) 

RROLL  = RROLL1  - RROLL2 

YYAW  =  TRNFRi0.4A  YW..\A3.BB3.KK22* RRL  DDE) 

SSP!  =  TRN FRO >.2. GGC30.CC3.DD3.EE i  ■ 

SSP2  =  rRNFRf  2.3.GGC40.CC4.DD4.KKC3  "EE2) 

SSP3  =  TRNFR(U,l.GGCiO.CC5.DD5.EE3i 
SSPEE1  =  TRNFR(  1.3.SSP.AA4.BB4.KK3 1  *FFIN) 

SSPEE2  =  TRNFR(0.2.SSP..\A5.BB5.KK32*RRU  DDE) 

PPOVVrER  =  SSP3-.SSP2-SSP1 

SSPEE3  =  TRNFRo).  1  ..SSP.AA<\BB6.KK53*PPOWER'; 

SSPEED  =  SSPEE3-SSPHE2-5SPEE  i 
CONTROL  FINT1M  =  30. 

*RINT  1., FFIN, RRL  ODER, RROLL. YYAW, SSPEED 
SAVE  (Sli  0.1. ROLL. RROLL. YAW. YYAW. SPEED, SSPEED 

GRAPHIG1  Sl.DE  =  TEK61S.PO  =  0..5)  TIMEl  LE  =  8..UN  =  SEC)  YAW(LO=-2.LI  =  1. 

SC  =  .5. NT  =  8}.YYAW(LO=  -2.LI  =  4,SC=  5.NI  =  8) 

LABELiGl )COMPARISON  BETWEEN  ROBERT  AND  F.M 
LABEL(G1  )FOR  YAW  OUTPUT 


a 


M  '  J.  *  «■  .  •  ' 

•T./ »  .*.*•■ 


W  /'.*  V  > V>  /  v  .  -•  *  -  --  V. >.  / \  * 

*.  S  .  T.  *-  %  *»  *.  A.  A  ■-  _  ■  *>  -  v  JV  -  •  M  ^ 


N  A  A  \  "• 


GRAPH(G2/S1,DE  =  TEK618.PO  =  0,.5)  T1ME(LE=  8., UN- SEC)  ROLL(LO=  -2.LI  =  1 
SC  =  .5,NI  =  8),RR0LL(L0  -  -2, LI  =  4, SC  =  .5,NI  =  8) 
LABEL(G2)COMPARISON  BETWEEN  ROBERT  AND  F.M 
LA3EL(G2)FOR  ROLL  OUTPUT 

GRAPH(G3,'S1,DE  =  TEK618,PO  =  0,.5)  TIME(LE=  8..UN  =  SEC)  SPEED(LO  =  -2, LI  = 
SC  =  .5, NI  =  8),SSPEED(  LO  =  -2.LI  =  4, SC  =  .5, NT  =  8) 
LABEL(G3)COMPARISON  BETWEEN  ROBERT  AND  F.M 
LABEL(G3)FOR  SPEED  OUTPUT 
END 
STOP 


SO 
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